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1 Introduction 



In recent years much attention has been focused on an issue of the phase transitions 
in multicomponent fluid mixtures, especially in ionic fluids. In spite of significant 
progress in this field, such systems are far from being completely understood. The 
investigation of complex models is of great importance in understanding the na- 
ture of critical and phase behavior of real ionic fluids which demonstrate both the 
charge and size asymmetry. The powerful tools for the study of multicomponent 
continuous systems are those based on the functional methods. In many cases the 
partition function of multicomponent models (see, e.g. [H) can be re-expressed as a 
functional integral after performing the Hubbard-Stratonovich transformation 0, 0] , 
a simple device proposed in the 50ies. Nearly at the same time another method, the 
method of collective variables (CVs), that allows in a explicit way to construct a 
functional representation for many-particle interacting systems was developed 0, 0] . 
The method, proposed initially in the 1950s for the description of the classical 

charged many particle systems and developed later for the needs of the phase transi- 
tion theory , was in fact one of the first successful attempts to attack the 
problems of statistical physics using the functional integral representation. Recently, 
the rigorous scalar field KSSHE (Kac-Siegert-Stratonovich-Hubbard-Edwards) the- 



ory [Hi, Il2| . which uses the Stratonovich-Hubbard transformation, was developed 
to describe the phase equilibria in simple and ionic fluids. As was shown 
both groups of theories are in close relation. 

In the CV representation of simple (one-component) fluids was reexamined 
from the point of view of statistical field theory. Our goal here is to derive the 
exact functional representation for the grand canonical partition function of a non- 
homogeneous multicomponent fluid. We reformulate the method of CV in real 
space and derive the CV action that depends on two sets of scalar fields - fields {p a } 
connected to the densities of the ath species and fields {u a } conjugate to {p a }- We 
study the correlations between these fields as well as their relations to the density 
correlations of the fluid. 

The CV method is based on: (i) the concept of collective coordinates being ap- 
propriate for the physics of the system considered (see, for instance, [l5j]) and (ii) 
the integral identity allowing to derive an exact functional representation for the 
configurational Boltzmann factor. Being applied to the continuous system the CV 
method uses the idea of the reference system (RS), one of the basic ideas of the 



liquid state theory [161 ] . The idea consists in the splitting of an interparticle inter- 
action potential in two parts: the potential of short-range repulsion which describes 
the mutual impenetrability of the particles and the potential describing mainly the 
behaviour at moderate and large distances. The equilibrium properties of the sys- 
tem interacting via the short-range repulsion are assumed to be known. Therefore, 
this system can be regarded as the "reference" system. The remainder of the in- 
teraction is described in the phase space of CVs (collective coordinates). The fluid 
of hard spheres is most frequently used as the RS in the liquid state theory since 
its thermodynamic and structural properties are well known. In this paper we de- 
rive the functional representation for the grand canonical partition function of a 
multicomponent fluid which includes both short-range and long-range interactions. 
The paper is organized as follows. In Section 2 we obtain the exact expression for 
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the functional of the grand partition function of a multicomponent non- homogeneous 
mixture. Section 3 is devoted to the study of the correlations of CVs fields and their 
relation to the density correlation functions of a multicomponent fluid. In Section 4 
we formulate the perturbation theory. The MF level of the theory is considered in 
detail. 

2 The functional representation of the grand par- 
tition function 

2.1 The model 

Let us consider a classical m-component system consisting of N particles among 

which there exist Ni particles of species 1, N 2 particles of species 2, ...and N m 
particles of species m.The potential energy of the system is assumed to be of the 
form 

1 m N m N 

a,/3 i^j a=l i=l 

where U a p(rf,Tj) denotes the interaction potential of two particles and the second 
term is the potential energy due to external forces. 

We present the pair interaction potential U a p(rf, if) as 

Uaf>tf, if) = vlptf, if) + W^lf, if), (2) 

where v®p(rf, if) is a potential of a short-range repulsion that can be chosen as an 
interaction between two hard spheres of respective diameters a a and a p. We call 
the m-component system with the interaction v®p(rf,r?) a reference system (RS). 
The thermodynamic and structural properties of the RS are assumed to be known. 
Wap( r T, r j) is some potential which can describe both repulsion (e.d. soft repulsion) 
and attractive interactions. In general,^/? (rf, if ) can be presented in the form 

W a p(vf, if ) = wSp(if, if) + w^(if, if), 

where uv?g(r",if) and w^Jrf,rj) are repulsive and attractive parts of the inter- 
action potential w a p(rf,Tj). Since w aj a(r",rf) are arbitrary in the core, i.e.for 

r ^ Qafi = {?a + °"/3)/2, we assume that the w a p(rf,Tj) have been regularized 
in such a way that their Fourier transforms w a p(k) are well-behaved functions of 
fcj and that w a p(0) are finite quantities. We denote by fl the domain of volume V 
occupied by particles. 

We present the potential energy UN 1 ...N m as follows 

U Nl ...N m = V§L.N m + \(p a \w a p\pp) + ty a \p a ) - N a v s a) (3) 

where 

P«W = E 5 ( r - r ") ( 4 ) 
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is the microscopic density of the ath species in a given configuration and fp a {r) is 
some external one-body potential acting on particles of species a. The following 
notations are introduced in ([3]): Vjyf Nm is the contribution from a m-component 
RS, z/f is the self-energy of the ath species 

= ^aa(O). (5) 

In (J3J) we have also introduced Dirac's brackets notations 

(Pa\w a p\pp) = / dr°drf p a (r")w al3 (r*, rf)p^(rf ), 



In the above formulas summation over repeated indices is meant. 

The system under consideration is at equilibrium in the grand canonical (GC) 
ensemble, (3 = 1/ksT is the inverse temperature {ks Boltzmann constant), p a is 
the chemical potential of the ath species. Then, the GC partition function can be 
written as 



RS 

Ni...N m 



N!>0 ± N 2 >0 * N m >0 

^{p a \w a p\pp) + (V a \p a ) 



(6) 



where V a {r) = u a + z/f — /3-0 a (r) is the local chemical potential of the ath species 
and v a = (3p a — 31nA a , A" 1 = (2Tim a P~ 1 /h 2 ) 1 ^ 2 is the inverse de Broglie ther- 
mal wavelength, (dr) = n^dr^, dTjsr a = dr" drf . . . dr^ is the element of the 
configurational space of N particles. 

For a given volume V, H[i/ a ] is a function of the temperature T and a log-convex 
functional of the local chemical potentials v a {r a ). 

2.2 The collective variables representation 

We introduce the collective variable p Q (r) which describes the field of the number 
particle density of the ath species. To this end we use the identity 



(7) 



exp Q (pHp)^ = j Vpb T \p- p] exp Q (p\w\p)^j . 

In (J7j) the functional "delta function" 8? [p] is defined as [13] 

5 F \p\ = jvu exp (z Hp)) , (8) 

Using (jSJ) we can present the Boltzmann factor which does not include the RS 
interaction in the form 

exp Q (pa\w a p\pyn = J VpVu exp \^{p a \w a p\pp) 

+i (uj a \ {p a - p a })) . (9) 
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Inserting equation in the definition flB]) of the GC partition function E[i/ a ] one 
obtains 



S KM] = J V P ex P (Pa\Wap\Pl3)j J[{pa, "a}] 



(10) 



where the Jacobian 



J[{p a , "a}] = J Vu exp (i{u a \p a )) E RS {{is a - iuj a }} (11) 

allows one for the passage from the microscopic variables r a (the Cartesian coor- 
dinates of particles) to the collective variables p a (r) (fields of the number density 
of particles). In (TTTi) Hr.s[{z7q, — i^ a }] = Srs[^*] is the GC partition function of a 
m-component RS 

N!>0 N 2 >0 N m >0 ° 

xexp(-/3V^ 1 s .. iVm )n^(r J ), (12) 



where Z*(r) = Y[ a Z a (r) = exp(i/*(r)). is the activity of the species a associated 
with the dimensionless local chemical potential z^(r) = V a (r) — iu a (r). It should be 
noted that J[{p a) v a }} does not depend on the pair interaction w a p{^i, ij) but only 
on the GC partition function of the RS 5rs[Z*] which is supposed to be known. 

Equation fflOl) can also easily be recast in the form of a standard statistical field 
theory, i.e. as 

S [{*VH = J VpVu exp {-H[{u a ,p a ,u a })) , (13) 
where the action Ti.[{u a , p a , u a }} of the CV field theory reads as 

W [{v a ,p a , uj a }] = ^ (p a \w a p\pp) - i (u a \p a ) - In S RS [{V a - iu a }\ . (14) 

Functional integrals which enter the above-mentioned formulas can be given a 
precise meaning in the case where the domain Q is a cube of side L ( V = L 3 ) with 
periodic boundary conditions. This means that we restrict ourselves to fields p a ( r ) 
and u) a (r) which can be written as Fourier series 

= eikr ' ( 15 ) 

keA 



and 

L 3 

keA 

where A = (2tt/L) Z 3 is the reciprocal cubic lattice. The reality of p a (and u a ) 
implies that, for k ^ p_k, a = Pita ( u -k,a — wjj, where the star means complex 
conjugation. Then, the normalized functional measure T>p (and T>uj) is defined as 
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vp = nn4^ ( i7 ) 

dp k , a dp_ kjCr = 2 dKp k>a d3p kjCr , k ^ . (18) 
Equation (117)) can be rewritten as 

v _ TT dp , Q -|-r d^Pg, Q d^p q , a 

where the sum in the r.h.s. runs over only the half A* of all the vectors of the 
reciprocal lattice A. We have for T>u, respectively 

Vuj = TT ^ TT d ^ q ' Q d ^ q '° ■ (20) 

Now let us present the action fTHj) for the isotropic interaction potential w a p(r) 
as follows 



w[{^,p„,w Q }] = - $ a/ 3(fc)Pk, a p-k,/3 - 1 y^ ^k, a pk,« 

k « k 

-lnS RS [{^a - iWa}], (21) 



Here CV p kjQ , describes the kth mode of number density fluctuations of the ath 
species. ^ a p{k) = yW Q/ 3(k), where w a p(k) is the Fourier transform of the interaction 
potential w a p{r). 

In order to obtain another equivalent representation of the action Ti.[{u a , p a , u> a }] 
we first distinguish the chemical potential v ° of the particle of the species a in the 
RS. To this end we put 

V a - iu a (r) =v° a - i^l(r) (22) 

and obtain 

\u a (r) = Au a + iu' a (r) 
with Au a = V a — i>®. As a result, the action (or the Hamiltonian) ( l2Ti) can be written 



as 



H[{v a ,p a ,u a }} = Az/ a p k=0 , Q + ^^^$ Q/3 (A;)pk,aP-k 



a,f3 k 

-i Yl ^k,aPk,a - In S RS - i^}] . (23) 

a k 

We have obtained the exact representations of the grand partition function of a 
multicomponent system (eqs. f|T3|) - f|T4|) . ff2T|) and (1231 ) in terms of CVs p a (r), which 
are number density fields of the ath species particles (or fluctuation modes of the 
ath species number density) We also stress that p a (r) and u a {r) are two real scalar 
fields and that eqs. ffTBl - fTHl) (as well as eqs. (|2"T1) and ([23])) are valid for repulsive, 
attractive as well as arbitrary pair interactions. 
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Besides the usual GC average (A)gc of a dynamic variable we introduce statis- 
tical field averages of the type 

(A [{p a ,uj a }]) cv = E y VpDu A [{p a ,uj a }} exp (-ft [{v a , p a ,uJ a }}) , 

(24) 

where ^4. [{p Q ,o; a }] is a functional of the two CV fields p a and u a . 

3 Correlation functions 
3.1 General relations 

Let us write some important relations. First, according to 2(], 21] we introduce the 
ordinary and truncated (or connected) density correlation functions 

GS.. Q JK}](l,...,n) = (n^(i) 

\l=l / GC 

1 5 n 2[K}] 



E[u a ] 5u ai (l) . ..5u an (n) ' 

Gai - a " [KIK '-' ?l] " ^ ai (l)...^ n (n) ■ (25) 

Our notation emphasizes the fact that the correlation functions (connected and not 
connected) are functionals of the local chemical potential ^(r) and functions of the 
coordinates (1,2, ... ,n) = (r 1; r 2 , . . . , r n ). For the sake of simplicity, we omit below 
the notations which indicate the functional de pen dence of the correlation functions of 
v a (r). In standard textbooks of liquid theory [l6( the n-particle correlation functions 
are more frequently defined as functional derivatives of H or log E with respect to 
the activities Z a = exp(i> a ) rather than with respect to the local chemical potentials 

p ai . ..«„(!, 2,..., n) I 5 n E ^ 



Z* ai (l)Z* 2 (2) . 5^(1)^(2) . ..5Z* n (n 

Pa,.. a „(l,2,...,n) 1 6 n \nE 



• (27) 

Z* ai (l)Z* a2 (2)...Z* n (n) E8Z* ai (l)5Z* m (2)...8Z* n (n) 1 ; 

We also define the partial distribution function g ai _ an (l, 2, . . . , n) and the partial 
correlation function h ai an (l, 2, . . . , n) [io| 

a (12 n) = ^°i- a "( 1 ' 2 '---' Tl ) (or) 

ya 1 ...a n K 1 -! ■ ■ ■ 1 11 /-,\ / n \ / \ l ZC V 

Pa 1 (l)Pa 2 (2)...Pa n (n) 

and 

h ai ... an (l, 2, . . . , n) = 5 1 i r ( ^ 2, """ , 1 V < 29 ) 

Expressions ( |25|) and fl28|) - (l29|) differ by the terms involving products of delta 
functions. For instance, for n = 2 and for a homogeneous system one has 

G$(l,2) = PaPl3 g a p(l,2) + p a 5 a p5(l,2), 

GQ' T (1,2) = Pa pph a p(l,2) + p a 5 al3 5(l,2), (30) 
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where p a is the equilibrium density of the species a. Besides, functions h a p(r) and 
9ap{ r ) are connected by the relation h a p{r) = g a p{r) — 1 for a homogeneous case. 



3.2 CV field correlations 



Let us consider the correlations of fields p a and u a . We start with the definitions of 
the correlation functions 



G£... P JR}](l,...,n) 
G£t..„JR}](l,...,n) 



n p«i (*) 

v«=l / CV 



n 



(31) 



cv 



and their truncated (connected) parts 



G p n) ' T p (1, 

Pa-, •■■f->a n \ 1 



G (n),T Q 



n 



n) 



m<n 



oJq,, ...u„„ V ' • " • ' 



n 



e n e£j 



(,^1) • • • ) ^mj ) 



(32) 



where the sum of products is carried out over all possible partitions of the set 
(1, . . . , n) into subsets of cardinality m < n. 

Correlation functions G p n } ... Pan . Now we introduce the modified partition func- 
tion 

S 1 [{v a , J a }] = J VpVu exp (-Hcv {{v a ,pa,u a }] + (J a \Pa)) , (33) 
where J a is a real scalar field and S 1 [{u a , J a }] is the generator of field correlation 



functions G p n J ... Po 



19j. As a result, we have 

1 5 n E l [{u a ,J a }] 



H^K}] 5J ai (l)...5J an (n) 
5 n logE 1 [{u a ,J a }] 



J Q „=o 



(34) 



The simplest way to obtain the relations between the G p } Pan [{z/ Q }](l, . . . ,n) and 
the density correlation functions is to start from the definition fT2o^ . One has 



Gg.. Q „(l,...,n) 



1 



5" S[{i/«}] 



H[{z/ a }] <Ji/ ai (l) . . .5fa„(n) 

5" H RS [{z7 Q - iu a }} 



+i (W a \p a )) 



5u ai {l) . ..8v an {n) 



sR^TI f vpVuj exp Q (p«\ w *<*fi\pp) 



x fir 



+i (u a \p a )) 

5 n E RS [{V a - iu a }] 
Su ai (l) . . .8u an {n) 



where we introduce the notation w* a[3 = —(3w a p and use the equality 
5 n 2 RS [R - iu a }] _ n 5 n S RS [{z7 a - iu a }] 



5u n , (I) ■■■ 5is n „ (n) 



Performing now n integral by parts yields 



G {n) (1 n) 

^ai...a n \ L } ■ ■ ■ , it>J 



S[R}] 



Su ai (l) . ..5u an {n)' 



VpVu exp ( - (p a \w* a/3 \pp) 



+ lnS RS [{z7 a - \uj a }}) 

I n 

n p<* (*) 

v»=l / CV 



5 n exp(iR|p a )) 
6u ai (l) . ..5u an (n) 



We have just proved the expected result 

G£U [K}] (1, • • • , n) = G^ pan [R}] (1, . . . , n) . 



(35) 



(36) 



Obviously, the following relation is also valid for the truncated (connected) correla- 
tion functions 



G$5* [R}] (1, - - - , n) = G^ pan [R}] (1, . . . , n) . 



(37) 



Correlation functions , Man . Let us define the modified partition function 

[R, J a }] = J VpVu exp (-Hcv [{v a ,p a ,u a }] + (J a \uj a )) , (38) 

where J a is a real scalar field. S 2 [{u a , J a }] is the generator of the functions Man 
and we thus have 



G { :1 .. Man [R}](l,...,n) 
GS^ an [{u a }](l,...,n) 



1 5 n Z 2 [{is a ,J a }} 



S 2 [R]} 5J ai (l)...5J an (n) 
5™log~ 2 [R,J Q }] 



SJ ai (l) . ..8J a in) 



(39) 



j a .=o 

In order to relate the correlation functions G$l Man and Ga™!..a n (l, ■ ■ ■ ,n) we per- 
form the change of variables p a — > Pq, + iJ Q in eq. (1381) . The functional Jacobian of 
the transformation is of course equal to unity and one obtains the relation 



1 



InS 2 [R, J a }] = -- {J a \w* a p\Jp) + lnS 1 [{i 



V a ,lW a/3 * 



(40) 



where the star * means space convolution and S 1 is defined in ( 1331) . The idea is to 
perform now n successive functional derivatives of both sides of eq. (j4"0l) with respect 
to J a . Since it follows from the expression (JMj) that 



5 n logH 1 ^,™* * Jp}} 



8J ai (l) . ..5J an (n) 



i n w* 



J a ,=0 



1,1 )...w* 



On«J V ' 



n, n 



one obtains 



(41) 



(Wai(l)>GV 

cS;L 2 [K}](i,2) 



c^L^ [K}](i,-..,n) 



cv 
1,1'} 



= i< Qi ,(l,l')(p ai ,(l')} ( 

= -w* lQ2 (1.2) -<«!,(!. 

x < 2Q2 ,(2,2')G( 2 j i f pQ2/ [{^}](l',2'), 

x Ggf.. pa; [K}](l',...,n'), n>3. 



4 The perturbation theory 



4.1 Mean-field theory 

Let us consider the functional of the GC partition function (1T31) with the action 
given by eq. (12T1) for the case of an isotropic interaction potential w a p{r). At the 



MF level one has 19 



Smf[{^}] = exp(-H{{is a , p a ,u a }]), 
where, for p a and uJ a , the action is stationary, i.e. 



5 H [{U a , p a , UJ a }} 



5p a 



5 H [{Va,Pa,UJ a }] 



5iO a 



. 



(42) 



(43) 



Replacing the CV action by its expression (12T|) in eq. (1431 leads to implicit equations 
for p a and uJ Q : 



iuJ a (l) = $^(1,2)^(2), 



(44) 



where pf s [{^ a — icJ Q }](i) denotes the ath species number density of the RS fluid 
with the chemical potentials {v a — iuj a }. For a homogeneous system ( 1441) can be 
rewritten in the form 



Pa = pT = Pa S [{*<*-iuJa}}, 
lU a = pp$ a p(0), 



(45) 
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It follows from the stationary condition (1431) that the MF density is given by 
p MF [K}] (1) = S\nEM F [{u a}] = ^ _ (i) ^ 



Su a (l 

and that the MF grand potential reads 



ft 



lnS MF [{v a }] = lnS RS [{z/ a - iu a }] + - (Pa F \w a p\pf F ) 



(47) 



The MF Kohn-Scham free energy of a multicomponent system defined as the Leg- 
endre transform 

PAmf [{pa}] = sup {(p a \u a ) - In E MF [{u a }]} , (48) 
has the following form in the MF approximation 

ftA M ¥ [{Pa}} = PArs [{Pa}} + ^ {Pa\w a p\pi}) - ~ I dr uw(0)p Q (r) . (49) 



Using the formulas 



G^ aia2 [R}] (1, 2) 
cfiL* [{P«}](1,2) 



£ 2 In S M f [{^a}] 
5z/ ai (l)5z/ Q2 (2) ' 

5 2 /3i M F [{Pa}] 



where (1, 2, . . . , n) = (ri, r 2 , . . . , r n ), one can get the well-known expressions for 
the partial pair correlation and vertex (or direct correlation) functions in the MF 
approximation. C^ F ai a 2 [{P"}] (1? 2) is obtained readily from the expression (1491) 



a 



(2) 



MF, apy^-i 



1,2) 



MF,o/3 l 1 , Z J 



where Crs^^I, 2) means the exact two-point proper vertex of the RS fluid at the 

mean field density p^ F . The two-point vertex function C^p a p * s connected to the 
usual direct correlation function of the theory of liquids cmf,o/3(1, 2) 



C MF,afl( 1 ' 2 ) - C M F,a/3(l,2) - 



Pa(l 



-<5 a /3<5(l, 2) 



In order to calculate Gy( F a Jl,2) we start with equation 



(1,2) 



df%*[{u a - iE7 }](l) <9p* s [R - iu7 Q }](l) 



0^(2) 



9^(2) 



(50) 



Pa S [{^«}] depends on directly but also through the mean field u a . Therefore, 
one has 



rr(2),T 



(1,2) 



^ s [K-iZU a }](l) 



X 



9^(2) 
<9u; 7 (3) 
9^(2). 



+ 



a P f[{^-iaJ a }](i; 



<9uJ 7 (3) 



(51) 
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Taking into account (|44p and (I50p we obtain finally 

Gmf^(1,2) = Gg^(l,2) -/3Gg } F T Q7 (l,3)^(3,4)G^(4,2). (52) 



(1521) can be rewritten in a matricial form as 12 



cgf (1,2) = Gif (1,2) -GS T (l,3M3,4)Gg T (4,2), (53) 

(2) T (21 r 

where Gj^/ R m(i,j) denotes the matrix of elements G^p,^ aj3 (i,j) and w(i,j) that 
of elements (3w a /3(i,j). The formal solution of eq. fl53|) is then 

GjJ F T (l,2) = (l + «;*^' r )~^G®' T (l,2) , (54) 
where 1 = 5 Q/ g5(l, 2) is the unit operator and the V denotes a convolution in space. 

4.2 Beyond the MF approximation 

In order to take into account fluctuations we present CVs p a and u a in the form: 

p a (l) = p a + Sp a (l), u a (l) =uJ a + Su a (l), 

where the quantities with a bar are given by (1441) . 

The function In E RS [{v a ] — i<^<*}] m CHD can be presented in the form of the 
cumulant expansion 

n>l ' ai,...,a n 

x9Jla 1 ...a n (l, • • • ,n)5u ai (i) . ..5u an (n), (55) 



where 9Jt ai ... an (l, . . . , n) is the nth cumulant defined by 

1 d n \nZ RS [{V a - iuj a }} 



9Jt (1 



'-i) n d5u ai {l) . ..d5u an (n) 



(56) 



8uj a ■ =0 



As is seen from (15 6 p and fl2"5"j) . the nth cumulant is equal to the n-particle partial 
truncated (connected) correlation function at u a = ZJ a . The expressions for the 
several cumulants given in the Cartesian coordinate phase space are as follows 

m ai (i) = Pai (i), (57) 

(l,2) + p Ql (l)<5 5(1,2)], (58) 
9Raia2ct3 (1) 2, 3) = p Ql (1)pq, 2 (2)p as (3)/i Q , lQ , 2Q , 3 (1, 2, 3) 

+p Ql (l)pa 2 (2)h ai a 2 (1, 2 )(5 5(1,3) 

+p Ql (l)p a3 (3)/i QlQ3 (l, 3)5 aiQ2 5(l, 2) 
+p Q2 (2)p a3 (3)/i 

0203 

(2,3)5 aiQ2 5(l,2) 
+p ai (l)5 aia2 5 aia3 5(l,2)5(l,3), (59) 
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"$^011012013014 (1> 2, 3, 4) p Ql (l)p a2 (2)pq, 3 (3)p a4 (4)/i Q , lQ , 2Q , 3Q , 4 (1, 2, 3, 4) 

i,j,k,l 

x5 am S(i,l) + P ai (^Pa k (k)h atak (i,k) 

i,j,k,l 

x5 aiaj 5 akai S(i,j)S(k,l) + Y Pai{i)p ai {l) 

i,j,k,l 

x h atCtl (i, l)5 ai0lj 5 aiOCk 5(i, j)5(i, k) + p ai (1) 
xtf ai£ *<y aia A ia4 <Kl, 2)5(1, 3)5(1, 4). (60) 

In the above formulas p ai (i) is the local density of the a«th species in the RS 
and h ai „. an (l, • • • , n.) is the n-particle partial correlation function of a m-component 
RS, defined in the GC ensemble (see eqs. fl^-®): 



h aia . 2 (l,2) = g aia2 (l,2) - g ai (l)g a2 (2), 

haxaiaz (1 , 2, 3) (y' Q , lQ , 2Q , 3 (1, 2, 3) 9a±a2 (1) 2)5 , a 3 (3) 

9 0.20.3 

(2,% ai (l) 

+2^(1)^(2)^(3) 
^"ctict2ct3c»4 (1, 2, 3, 4) ^q, iQ , 2Q , 3Q , 4 (1, 2, 3, 4) .... (61) 

In the case of a homogeneous system a Fourier image of the nth cumulant can be 
presented in the form 

SWax^Cki,...,^) = ((iV^),...^)) 1 ^,..^^,...,^) 

x5 kl+ ... +k „, (62) 

where 5' ai ... Qn (ki, . . . , k„) is the n-particle partial structure factor of the RS. 
Substituting ( l55|) in ( fl~3l) one can obtain 

= Srs [{ v a - J V5pV5u exp j-~ (5p Q |w Q/3 |(5p^) 
+i(^ Q |5p Q ) + ^^ ^ fdl...Jdn 

n>2 ai,...,a n 

xTl ai ... an (l,...,n)5uj ai (l)...5uj an ( y n)} . (63) 
Integrating in (1631 over <5co> a . (i) we have in the homogeneous case 

"KM] = 2 M f2' J (dSp) eX P { ~ ^| 5^ L «/3( /E )5Pk,a5P-k,/3 

+ ^ft„(5p Q )|. (64) 



n>3 
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Gaussian approximation. In the Gaussian approximation, which corresponds to 
taking into account in (1641 only the terms with n < 2 (H n = 0), we have L a p(k) = 
C a p{k), where C a p{k) are the Fourier transforms of the partial direct correlation 
functions. After integrating in ( |64l) we arrive at the GPF of a m-component system 
in the random phase approximation (RPA). 

Using the Gaussian averages one can develop a loop expansion of SK^oJ] in the 
CV representation as it was done recently for a one-component fluid [131 ]. 

5 Conclusion 

Using the CV method we have reconsidered the basic relations of statistical field 
theory of a multicomponent non-homogeneous fluids that follow from this approach. 
In contrary to the KSSHE theory [12] the corresponding CV action depends on 
two sets scalar fields - field p a connected to the number density of the ath species 
particles and field u a conjugate to p a . We derive the explicit expressions for the CV 
field correlations and obtain their relation to the density correlation functions of a 
multicomponent system. 

Contrary to the theories based on the Stratonovich-Hubbard transformation 0, 
0] , the CV representation has some important advantages which could be very useful 
for more complicate models of fluids. In particular, it is valid for an arbitrary pair 
potential (including a pair interaction w; a( g(l,2) which does not possess an inverse) 
and is easily generalized for the case of n-body interparticle interactions with n > 2. 
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